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In two different phase transition models of nuclear frag- 
mentation we show that the emission of pre-equilibrium par- 
ticles and mixing of events from different classes cannot be 
ignored in the analysis of nuclear fragmentation data in terms 
of critical exponents, and we show how the apparent values 
of the extracted exponents are affected. 



The premier goal of medium and high energy heavy-ion 
reactions is the exploration of the nuclear phase diagram. 
On theoretical grounds we expect infinite nuclear matter 
to undergo at least two distinct phase transitions. One 
is the deconfinement or quark-gluon-plasma phase tran- 
sition. The other is a "liquid-gas" type phase transition. 
It is believed to be of first order, terminating at the crit- 
ical point in a second order transition. In nuclear mul- 
tifragmentation reactions one attempts to map out the 
liquid-gas coexistence region and locate the critical point. 

The first dataset to be interpreted in terms of crit- 
ical exponents resulted from proton-induced spectator- 
fragmentation of krypton and xenon targets . How- 
ever, it was later shown in the framework of the per- 
colation model H that value of the "critical" exponent 
r observed in H0] was predominantly a result of mix- 
ing of different event classes and integration over impact 
parameter. In addition, similar power-law behavior was 
seen in classical reaction dynamics simulations, where it 
was shown to be inconsistent with matter going through 
the critical point ||-|(| . 

A significant step forward was then taken by perform- 
ing event-by-event analysis of the moments of the mass- 
or charge-distributions of the fragments • The result 
of this analysis of emulsion data suggested that nuclei 
break up similar to percolation clusters |t]]. 

It was hoped that the problem of impact parameter 
selection would be less severe and therefore the critical 
exponents easier to extract for participant fragmentation 
in symmetric heavy ion collisions. By focussing on very 
central collisions as a function of beam energy, a mini- 
mum of the fit parameter A, with cr(Zf) cx ZJ X was ob- 
served |J. This confirmed similar observations obtained 
from reverse-kinematics reactions p~QQ . It had been pre- 
dicted that this minimum would correspond to the actual 
value of the critical exponent r gj. However, theoretical 
calculations showed that there may be significant topol- 
ogy effects such as the formation of bubbles and toroids 
|pr| |l3| at work, changing the values of the measured ex- 
ponents [Q. 



At present, significant effort is also directed at the 
study of the moments of the fragment charge distribu- 
tion in reverse-kinematics reactions, where the rise and 
fall of multifragmentation was observed (l^Jl^,|l^]. Re- 
cently the EOS TPC collaboration used the reaction 1 
A GeV Au + C in an attempt to reach the critical point 
of nuclear matter and determine the critical exponents 
in the spectator fragmentation of the residue of the gold 
nucleus 

It is generally agreed that the fragmentation process 
in proton-induced or reverse-kinematics reactions pro- 
ceeds in two steps - a first pre-equilibrium step in which 
the participants interact and deposit excitation energy 
into the spectators, and a second equilibrium step in 
which the excited (and hopefully equilibrated!) specta- 
tor residue decays. To this end we have constructed a 
hybrid model, in which the pre-equilibrium energy de- 
position and resulting residue size are calculated in the 
framework of an intra-nuclear cascade (INC) model Jl9| , 
and in which we calculate the decay of the residue within 
a percolation model H jq| or statistical multifragmenta- 
tion model (SMM) |2(J framework. For any given im- 
pact parameter the INC provides the charge and mass as 
well as the excitation energy per nucleon of the spectator 
residue by calculating a sequence of individual nucleon- 
nucleon collisions and single-particle removals from the 
nuclear potential well. It is worth noting that in cen- 
tral collisions the average residue charge is less than 49, 
implying the emission of 30 pre-equilibrium charges. 

The excitation energy can be converted into a perco- 
lation bond-breaking probability via 0] 
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where T is the generalized incomplete gamma function, 
B is the binding energy per nucleon in the residue (taken 
as 6 MeV here), and T is the temperature as calculated 
from the INC model, T = y/E*/a. The SMM is micro- 
canonical statistical multifragmentation model where the 
decay probabilities into channel j is proportional to the 
statistical weight Wj cx exp[Sj(T, A, Z)], where Sj is the 
entropy in channel j. In the SMM we used a standard 
set of parameters applied before for analysis of multi- 
fragmentation in peripheral heavy-ion collisions fig] . In 
particular the break-up density is assumed to be 1/3 of 
the normal nuclear density. 

We then arrive at the final inclusive mass distribution 
(and its moments) as follows: For any given event we 
select an impact parameter with geometrical probabil- 
ity weighting, calculate E* — > T — > pb and the pre- 
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equilibrium multiplicity, m pro , within the INC stage, and 
then use this ph and residue charge, Z ICS — 79 — rn prc , 
to obtain the multiplicity, m eq , (and charge/mass dis- 
tribution) of the equilibrium particles from the percola- 
tion stage or SMM stage. The total multiplicity is then 
A standard integration over impact 



fc prc 



parameter yields the desired distributions. Within the 
assumptions stated above, both hybrid models are then 
free of adjustable parameters. 

In Fig. 1 we compare the results of our models to the 
data of |l8| . Displayed is the second moment of the 
charge distribution, computed event-by-event, and aver- 
aged over all events with identical total charged particle 
multiplicity, 



M 2 {m) = y^<5 m , mi 
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Here Ni{Z) is the number of fragments of charge Z emit- 
ted in event i and rrii is the total charged particle multi- 
plicity of this event. is a summation over all events 
and S m . mi is the usual Kronecker-<5. Z max is the upper 
cutoff used in the summation over all fragment charges. 
Here we should point out the usefulness of displaying 
multiplicity-sorted data of M 2 (m) by including and ex- 
cluding the largest, 2-largest, 3 rd -largest, . . . fragment. 
This information is available in a model-independent way, 
and it shows the most probable partition of the system 
in each multiplicity bin. In the EOS data |Q displayed 
here (circles) only part of this important information is 
given: The largest fragment was included for multiplic- 
ities above 26 and excluded below. For the models we 
display in the lower curves M2(m) excluding the largest 
fragment and in the higher curves including the largest 
fragment for all multiplicities. Fission events were elimi- 
nated from both the data and the calculations. 

We observe an astonishing degree of agreement be- 
tween the results of the INC/percolation hybrid (his- 
togram) and the data for all multiplicities. We should 
point out here that the fragmentation of a Z = 79 residue 
with the same excitation energy distribution results in 
average moments very different from the data, indicat- 
ing the importance of the emission of pre-equilibrium 
charged particles for this observable. In addition we per- 
formed a calculation fragmenting a Z — 79 residue at 
a fixed excitation energy corresponding to the critical 
point of the percolation model. This resulted in events 
distributed over the multiplicity interval from 12 to 28, 
having approximately constant values of M2 ~ 100 when 
excluding and M2 ~ 3000 when including the largest 
fragment. This also indicates the importance of mixing 
between events of different excitation energy classes in 
the observed values of the moments. 

The INC/SMM hybrid shows also qualitatively the 
same features as the data, but is generally too high on 
the lower branch and too low on the upper branch by 



about a factor of two. Presumably one could generate 
better agreement for the INC/SMM by fitting an exci- 
tation energy and residue size distribution to the data 
along the method used in ]l6]] . This, however, is not the 
goal of the present work. 

Rather we would like to point out that the data and 
both models have similar slopes of M 2 versus m over the 
entire multiplicity interval. These slopes were used by 
the EOS collaboration to extract the "critical" exponent 
7. This is even more obvious in Fig. 2, where we display 
logAf2(m) as a function of \m — m c \ with m c = 26, ac- 
cording to the analysis of |l7]] . We see that data and both 
models have roughly parallel "liquid" and "gas" branches 
over the interval displayed here. However, in contradic- 
tion to the analysis of |^l[ , we do not recover the known 
critical exponent 7 (=1.8) of the percolation model from 
this figure. Instead both models yield about '7=1.4, in 
accordance with the experimental observation. 

We have also performed an analysis of the size of 
largest fragment as a function of the total multiplicity. 
From an analysis like this, the authors of ref. |17| claim 
to have extracted the critical exponent (3. Recently, one 
of us |22| pointed out that the presence of pre-equilibrium 
particles may change the value of the extracted apparent 
exponent. Within the framework of the present hybrid 
models we can now perform a quantitative comparison 
to the data, as shown in Fig. 3. Performing a linear re- 
gression fit to the output of the INC /percolation model 
in the interval shown yields a slope of '/3'=0.35, reason- 
ably close to the value of 0.29±0.02 extracted from the 
data, and not close to the nominal value of /3=0.41 for 
the percolation model. A fit to the INC/SMM output 
yields '/3'=0.50. 

It is dangerous to interpret these exponents extracted 
in the above way as the true critical exponents charac- 
terizing the universality class of the actual nuclear phase 
transition. Three effects play important roles in this: the 
extreme finite size of the fragmenting system, the use of 
log \ m — m c \ instead of log \ T — T c \ (or log \p — p c \) for 
the abscissa, and, most importantly, the mixing of dif- 
ferent event classes of size and excitation in the same 
multiplicity bin. This mixing is caused by the use of in- 
clusive, impact parameter integrated data with different 
admixtures of pre-equilibrium particles. 

The mixture of events with different "temperature" 
into the same multiplicity bins represents a convolution 
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where 0(m) and 0'(T) is the physical observable as a 
function of multiplicity and temperature, and N(m,T) 
is the (integral normalized) number distribution func- 
tion of events with a given temperature and multiplic- 
ity. In Fig. 4, we show N(m, T) as a contour plot for the 
INC/percolation model. Each contour level represents 
a factor of 2 higher value than the one surrounding it. 
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One can clearly see that N(m, T) is rather broad, indi- 
cating a large degree of mixing of events with different 
temperature into the same multiplicity bin. 

While the scaling laws of critical behavior are valid for 
0'(T) in the vicinity of the critical point, no such state- 
ment can a priori be made for the same observable as 
a function of the multiplicity. To accomplish this, the 
inverse of the convolution function, N(m, T) _1 has to be 
known. For obvious reasons, the function N(m,T) can- 
not be extracted from data in a model-independent way, 
and therefore the inverse is also not accessible. For all 
these reasons, a precise model-independent determina- 
tion of critical exponents, as sugg estedinrefs. 
is not correct. 

In conclusion, taking the data of the EOS collabora- 
tion we have shown that an interpretation in terms of 
critical phenomena should be carried on very carefully. 
Mixing of different event classes into the same multiplic- 
ity bins and the contributions of pre-equilibrium particle 
emission make a mo del- independent extraction of critical 
exponents questionable. 

What is possible, however, is a detailed comparison of 
the data to models which include pre-equilibrium as well 
as equilibrium components, and which emulate as closely 
as possible the experimental trigger conditions. We have 
used two of the models here to perform this kind of com- 
parison. We have found, e.g., that the INC/percolation 
model achieves surprising agreement with the data. An 
analysis of the kind employed in for the model out- 
put yields apparent exponents of '/3'=0.35 and l 7'=1.4, 
whereas the known critical exponents for 3d-percolation 
in the infinite size limit are /3=0.41 and 7=1.8. The 
next step, in our opinion, could be analysis of exclusive 
events aimed at finding actual equilibrated residues. Af- 
ter that we can really investigate manifestations of the 
phase transition in the residue-nuclei. There is hope that 
this direction of research will be fruitful. The rather sur- 
prising agreement of the INC/percolation model (which 
explicitly contains a phase transition) with the data given 
rise to this hope. 
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FIGURE CAPTIONS 

Fig. 1 Second moment of the charge distribution as 
a function of the charged particle multiplicity 
for the reaction 1 AGcV Au + C. Histogram: 
INC/percolation, line: INC/SMM, circles: data 




Fig. 2 Second moment of the charge distribution as 
a function of the charged particle multiplicity 
for the reaction 1 AGeV Au + C. Histogram: 
INC/percolation, line: INC/SMM, circles: data 




Fig. 3 Charge of the largest fragment as a function of 
the charged-particle multiplicity for the reaction 1 
AGeV Au + C. Histogram: INC/percolation, line: 
INC/SMM, circles: data 0. 

Fig. 4 Histogram of the number of events as a func- 
tion of the temperature, T, and the total charged 
particle multiplicity, to, as predicted by the 
INC/percolation model. Each contour line is off- 
set by a factor of 2 from the previous one. 
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